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ABSTRACT 


The algorithm presented is an extension of the Land and 
Doig branch and bound method combined with the’ branch 
selection techniques presented by Beale and Small to solve 
integer or mixed-integer linear programs. The algorithm 
obtains the solution by solving a linear program with upper 
and/or lower bounds on selected branch variables. By 
systematically changing these bounds, and maintaining only 
the current canonical form, the solution is assured using a 
minimum of excess computer storage above that required to 
solve the linear programming problem. Thus the problem can 
be solved entirely within the computer core, and the problem 
converges to the solution faster than most other’ general 


integer linear programming algorithms. 
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1, INTRODUCTION 


In the literature today, there are many” algorithms to 
solve Integer or mixed-integer linear programs, but most of 
these algorithms are for spectalized types of problems, such 
as transportation problems or machine scheduling problems. 
As Bellmore and Nemhauser [2] say, "Integer linear 
programming algorithms are notorftous for converging rapidly 
On one problem and then performing miserably on the next.” 

The algorithm presented here fs Harold Greenberg's 
natural extension of the Land and Doig branch and hound 
method [s] ; [9] combined with branch selection techniques 
presented by E.M.L. Beale and R.E. Small pu] , to solve 
Integer Linear Programs. In additton, the algorithm has 
heen developed to solve mixed-integer linear programs where 
some, but not all of the variables are integer-restricted. 

The problem is solved hy primal linear programming 
using upper and lower hounds on an intelligently selected 
integer-restricted variable that does not have an Integer 
value In the current solution. These hounds are’ system- 
atically changed until an optimal solution [fs assured. In 
addition to speeding convergence by having less’7 branches 
than most other methods, this method also uses a minimum of 
excess computer storage ahove that required to solve’ the 
linear programming problem. Thus the need for offline 
storage is eliminated and the enttre problem fs solved and 


maintained In the core, giving a relatively small cpu time. 


Computational results and comparisons with other algorithms 


are Included to show this small time. 


ij, FORMULATION OF THE PROBLEM 


The problem to be solved is: find X = (x4,...,x,) that 
minimizes  Cx' 


(1) 
when AX'=B8 


0< Xj G Uj 2) f=l, oe ia 
a integer for some or all jzl,...,n. 
C is an n-vector, B is an me-evector, and A is an mxn matrix, 
and the 2 are positive integers (or are infinite). 
To formulate a problem into this form, one need not 
restrict the slack/surplus variables to being integer, thus 
there is no need to restrict the elements of the A, B, or C 


vectors to being integer. 


The solution to (1) includes the case where Xj is 
restricted to being 0 or i. This accomplished by 
restricting ay to being integer with an upper bound of one 


(i.e. the corresponding a equals one). 


TH RA ING PROCEDUR 
A. THE LAND AND DOIG METHOD 


The branching procedure employed In this: algorithm Is 
an extension of the branch and bound method of Land and Doig 
191; their method will be explained first. Land and Dolg 
first solve (1) as a lInear program to obtain the contInuous 
solution, (2° > Mi Piers «90 ae If the continuous’ solution Its 
all Integer, the problem Is solved. tf there are fractions, 
this iInitlal value of the objective function, 2, serves as 
a lower bound on_ the value of any feasible optimal 
Integer-solution. One variable, Xj Is then selected whose 
value is not Integer. In the final Integer-solution, Xj» 
must satlisify elther Xj <5 Jum 0 Xj > (xj ) + 1, where 
(x5 ) ts the greatest Integer less than Xp The procedure 
initlally Involves the soluttion of two problems: both 
contain the Initial linear programs, but with an additional 
constraint; In one case the constratnt Is SUNS (x7 J. and 
in the other case it Is Xj = x7 J} + 1. Once these 
solutions are obtained, and supposing that both are _ not 
integer-solutlions, the one with the minimum objective value 
Is Investigated first, while the other Is stored for later 


Investigatlon. Suppose the minimum of the two objective 


functions Is zi , With solution: 
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(z 1 


eee, (x Ob... ,x*). Now a new vartlable, x, , whose 
1 J n k 

value is not integer, Is selected, and three new problems 

are solved: these are the initial l!near program (1) subject 


to the additional constraints x, = (x; J, X, * Bes J, tn the 


J 
first case, X j = [x9 ale Xs xi J + 1, In the second 
case, and x, = cx? J - 1, In the third case, all assuming 


[x J #0. 

The procedure is again repeated until an integer- 
solution ts obtained. At every step where there is a new 
varlable chosen for restriction, there fs an additional 
problem to solve. Thus the next step would have four 
problems to solve, and the one after that, five, If the 
Integer-solution ts not obtained. Once an integer-soluttion 
is achieved, its objective function value, 24 , becomes a 
new upper bound on the optimal solution. With this bound, 
the procedure then Investigates the other problems that are 
solved at this step. At any point in the’ branching, when 
the current upper hound is exceeded, or when an 
integer-solution is obtained, that branch is terminated. 

This complete process enumerates the possibilittes for 
the solution to (1) ina directed tree. The rooted node 
corresponds to 2° and ft has directed branches to the nodes 
corresponding to the solutton of the Jinear program (1), 
with x5 set at a specific value, each of these’ integer 
values representing a single’ branch. At each of these 
nodes, the branching ts repeated. The mintmum objective 
value at each node determines the next branch, which limits 


the potentially large size of the tree. 
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The problems that are stored for further Investigation 
rapidly create an Immense storage prohlem, thus the Land and 
Dotg method becomes unwieldy for problems with more’~ than 
Just a few variables. 

The algorithm presented here takes the Initial thought 
of the Land and Doig method, but uses an_ inequality 
constraint for the selected variable. This inequality may 
be aS an upper or a lower bound for the variable. By 
employing a bounded variable linear programming solution 
technique similar to that from Dantzig {3 ), the algorithm 
allows the problem to be solved with a minimum of excess 
storage used above that storage required to solve the linear 


programming problem. 


B. THE NEW BRANCHING PROCEDURE 


As with the Land and Doig method, the procedure obtains 
the continuous solution to the problem. If this solution is 
not an Integer solution, the procedure then selects a 
non-integer valued varlable, Xj , to branch upon. This 
creates two new problems; both contain the original linear 
program (1), but one has’ the added constraint of xj OX] 
while the other has the added constraint of i Cx") a 
Inittally, only one of these problems fs_ solved. The 
solution to this problem is treated In the same manner, and 
the process Is continued until] an Integer solution is 


obtained, or the problem becomes infeasible. One thing to 


notice Is at every node there are only two. possible 
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branches, and there are no problems saved. This limiting of 
branches, alone, tends to reduce the number of nodes 
required for the solutton of the problem. In addition, once 
an Integer solution Is obtained, it serves as an upper bound 
on the final solution, and no branching continues beyond a 
node where the objective value excedes the current upper 
bound solution, 

The strategy used ts to Investigate a single branch of 


the tree until z°? 


> zz , where zz ts the current upper 
bound on the optimal solution, At this point, backtracking 
begins, to find a node that has only one previous branching. 
This backtracking ts accomplished simply by removing. all 
bounds from the nodes that are further out on the tree, and 
solving the resulting problem. This solutton ts the 
re-creation of a previous solution that would have’ been 
stored otherwise. The second branch ts then made utllizing 
the Initial problem at that node which Is re-created from 
the latest transformation of the equations, and the 
constraint which was Initially tgnored at that node, and 
solving this linear program. At times, the procedure may 
impose further bounds on aé_e previously bounded vartable. 
This presents no problem since’ the inittal bounds will be 
satisified. 

Thus, this procedure stores only the current solution 
and the current bounds. The optimal solution ts_ the 
existing upper bound solution when’ It becomes Impossible to 
branch any further, This ts realized when’ backtracking to 


the rooted node a second time during the solution procedure, 
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V 0 UR US 


The program utIlized the two phase simplex method to 
solve (1), together with a bounded varlable’ technique 
similar to that In (3). By using this technique, the 
problem matntains m constraint equations throughout’ the 
solution procedure and has a_ solution where’ some of the 
varlables willl be at thelr upper bounds, some of the 
varfables will be bastic, and some of the varlables will be 
at thelr lower bounds. 

At each node In the procedure, we are Interested In the 
solution to (1) which satisifies the new’ bounds given by 


O41, Coxe vie < 
J J J _ 
upper bound, which may be Infinite, uj Is the Imposed upper 


u., jzl,...,n, where u. Is the problem 
J 


bound generated by the branching procedure, and 1, is the 
Imposed lower bound, also generated by the branching 
procedure. After achleving a solution for the “s which 
satisifles the constraints and the Imposed bounds, we can 
obtain the following canonical form [3], from the constraint 


equations in (1): 


n 
X, +8 2 Yonex. =d., t=1,...,m, 
jenn ie ee 
(3) 
n 
=2 Geo C ee. = -2, 
j=m+i1 J J 


where the x. , t#l,...,m, are the basic varltables, and the 
x , J=m+l,...,n, are the non-basic variables. When there Is 
J 


a feasible solution to this problem (3), with upper = and 
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lower bounds, we can use a_ varlation of the usual simplex 
method to obtain optimality. This variation is suggested hy 
utllizing a method to tmprove a feasible solution’ that is 
not optimal. To tmprove this feastble solution, we can 
increase a non-basic variable x, at its lower bound when 
© <0, or we can decrease a non-hasic variable x at its 
upper hound when - >» 0, where 2 is the or. cost 
coefficient for the ; th variable as found In the ohjective 
function In the cononical form. These conditions’ for 
optimality come from the following theorem |5]|: 
Theorem: 

any values of X jzl,...,n, satisifying the constraint 


equations in (3) and 0<¢ 17. < x | TT Oe 


ol li 7 pelt 
are an optimal solution for minimum z, with objective 


vaiue z , Iff Cc. > 0 for variables at thelr lower bound, 
O J 


T., and c. < 0 for vartables at thelr upper bound, u 


j j 


The proof is obvious, for any tincrease tn a variable at tts 


lower bound, or any decrease in avartiable at Its upper 
bound an only Increase z. More detailed proofs may be 


found tr |3], or In [6]. 
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V, THE VARIABLE SELECTION PROCEDURE 


While any non-integer varlable [In the present solution 
can be chosen for the branching, the size of the directed 
tree can be further reduced by selecting the one which will 
produce the optimal solutton quickest, thus reducing the 
number of branchings and nodes. The method employed [In this 
algorithm was developed by Beale and Small [4], using a 
suggestion of Driebeck. 

The process Involves comparisons of the marginal cost of 
Increasing and of decreasing each of the non-Iinteger basic 
variables In the current solution, as in the cost ranging 
procedure of linear programming, where the marginal cost of 
changing each of the non-basic varlables Is considered, and 
the non-hasic vartable with the greatest change which still 
allow feasibility Is entered Into the basis. The marginal 
cost, D,., of Increasing X ; is given by 

D 5; hey . Cc, Aer 
VJ 
and the marginal cost, D)., of decreasing x, Is given by 


. = Min (c. / y.. ). 
Do yy... > 0 S “i 
Vy 


where C; , Is the reduced cost coefictent for the 


Le 


jth 
variable as found In the objective function In the canonical 
form. By applying the strategy of taking the varlable that 
makes the worst of the two alternatives as bad as possible, 
we have the alternative that Is the most unlikely to contain 


the final solution, and we branch on this variable in the 
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other direction. 

This strategy is tmplimented by computingd ,. D,. x fi 
and ey; aD .* (1 - f,), for each non-integer basic 
varlable that Is restricted to being Integer. Where f. is 
the fractional part of xX, 4 X.5 [x.] + f.. The branching its 
then made on the varlable that makes the larger of these two 
quantities as large as possible. Say that Max J. > 
Max yj. we then go hack to find i such thatL, = es, ° 
Thus, the pth bastc vartable would be the one branched upon, 
and the direction of the branching ts x. 2 {x, ] + 1, since 


It Is the/. Dj and not the. U3 that Is the greater. 
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V T OR'THM 


A. THE BRANCHING PROCEDURE 


The branching procedure ts contaltned In the following 


algorithm: 
1, Definitions 


S(ZZ,X) oeaeeX) is the current upper bound solution 
which represents a_ feasible solution to (1) with objective 
value zz, and having all varlables at Integer values which 


are required to be Integer In the solution. 


h bounded vartable. 


B(L) Is the bound restricting the pth bounded 


i(L) Is the Index of the i* 


variable, 

BSIGN(L) Indicates whether the bound on the he’ 
bounded variable Is an upper or lower bound: BSIGN(L) = + 1 
for an upper bound, and BSIGN(L) = - 1, for a lower bound. 


N(L) Is the counter to show the number of branchtiIngs 


for the tt hounded varlable; Initially, NCL) = 0. 
2 Logic Proce 


a. Set L=0, and take zz as’ Infinite, unless there 
Is some known feasible Integer-solution to (1), to get the 
Initial S(2z,Xx, peeeeX). Solve (1) as a linear’ program to 
obtain the continuous’ solution. If this solution has 
Integer values for all of those variables that are Integer 


restricted, the problem ts_ solved. Otherwise go to ob, 
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maintaining the canonical form of the solution to (1). 
b. Set L=L + 1. Go to b(1). 

(1) Select the hasic variable, x , that is 
fractional, but required to be integer, which is the most 
unlikely to produce the optimal solution, and go to b(2). 

(2) Set BSIGN(L) opposite the way which will 
produce the most unlikely chance for the optimal solution 


and go to b(3) or b(&). 


(3) If BSIGN(L) = +1, set B(L) = - . Go to 
| om 

(4) If BSIGN(L) = - 1, Set B(L) = - + 1. Go 
wO Cc. 


c. Solve the linear program using the maintained 
canonical form with the new bound on the variable 1G 
where j = I(L). Go to the following cases, one of which 
will hold: 

(1) !'f the solution produces an objective 
value, z, with z > zz, go to d. 

(2) If the solution produces an objective value 
z, with z < zz, and the solution has integers for all of the 
required variables, redefine S(ZZpX poe need ) as a new 


n 


feasible tnteger solution upper bound, where 2z = Z, Ae 


Bein? is the new solution. Go to d. 

(3) If the solution produces an objective 
value z, with z < 2z, and the solution does not have 
integers for all of the integer-restricted variables, go to 


b. 


a, 


(4) If the problem has an InfeasIible solution, 
£0 to d, 

d. Set LL = L, and to to d(1l). 

(1) If N(L)=2, go to d(2). Otherwise N(L)=1; 
POutO dS). 

(2) If L=l, the current feasible solutton 
eo / gee Jis optimal, and_ the algorithm stops. 
Otherwise, set L=L - 1 and go to d(1l). 

(3) Solve the linear programming problem 
Starting from the current canonical form with the bounds 
BCL), BCL + 1),...,B(LL), removed. tf j = I(L), and x; does 
not become a basic variable, the solution Is non-unlque; sur 
is then made a basic varlable. Go to d(4&) or d(5). 

(4) If BSIGN(L) = +1, set BCL) = BCL) +1, 
BSIGN(L) = <1, and N(L) = 2, and go to c. 

(5) If BSIGN(L) = -1, set BCL) = BCL) ol, 
BSIGN(L) = #1, and N(L) = 2, then go toc. 


This completes the branching part of the algorithm. 


B. THE VARIABLE SELECTION PROCEDURE 
The selectton procedure for which basic variable’ to 
branch upon, and which way to branch Is contalned In the 


following algorithm. 


1, Definitions 


! ts the ‘aa basic varlable. 


h 
J is the ae column In the basis. 


IB(1l) ts which problem vartable Is_ the jth baslic 


varlable. 
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IMU Is the Index for which |! has the maxd. yj : 

IML Is the Index for which | has the maxd yp, 

DUC(J) Is the current ratio = lc, /y, Lot vi S Or, 

DL(J) Is the current ratio = loj/y, jl iF vin 0. 

XDU Is the current value of the up DUCJ). 

XDL Is the current value of the min DL (J). 

XLDU Is the current value of the max ia 

XLDL Is the current value of the max nie 

BSIGN Is the f[ndicator for whether’ the’ branch 
requires an upper bound or a_ lower bound on the varlable to 
he restricted. BSIGN = +1, for an upper bound, and BSIGN = 
-1, for a lower bound. 

M is the number of constraint equations, tf.e. the 


size of the basIs. 


2, logic Procedure 


a. Set XLDU=XLDL=IMU=IML= 0, I= 0, and go to b. 
b. Set f=! + I. If IBCl) Is an artificial 


variable, or If Is not restricted to belng Integer, 


*1B(1) 


or If is Integer, go back to a. Otherwise, set J=0, 


*1BC1) 
XDU=XDL= Infinity, and go to c. 


c. Set J=JS+1, If = = 0, repeat step cae =If 


v3S 0, go to c(1), If Yay > 0, go to c(2), otherwise, Yh 
0, and repeat step c. 

(1) DUC) = Ic ;/y 41. If DU(J) < XDU, set XDU 
= DU(J), and go to c(3). Otherwise, go stralght to c(3), 


without setting these values. 


(2) DLJ) = le ;/yj jl. If DL(J) < XDL, set XDL 
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= DL(J), and go to c(3). Otherwise, go straight to c(3), 
without setting these values, 
(3) If J<M, go toc, otherwise this is the 
last comparison for the 1th basic variable and go to d. 
d. Go to the following cases. 
(1) «f XDL is still equal to infinity, XLDOL its 
unchanged, otherwise, if XDL x f >» XLDL, set XLDL = 


IBCI) 


XDilax Figct) and set IML = I, go to d(2). 
(2) tf XDU is still equal to infinity, XLDU is 


unchanged, otherwise, if XDU x (1 =- Fipcry? > XLDU, set XLDU- 


= XDU x (1 ), and set IMU = 1, go to d(3). 


“FB 
(3) If 1t> MM, go to b, otherwise this its the 
last basic variable to consider, and go to e. 
e. One of the following cases will hold. 

(1) If IMU and IML both equal zero, this 
procedure was called in error, and there are no variables to 
branch upon. 

(2) If XLDU > XLDL, then the variable to branch 
upon is IMU, and BSIGN = + 1, to set an upper bound on Xemu’ 

(3) «tf XLDU < XLDL, then the variable to branch 
upon is IML, and BSIGN = - 1, to set an lower bound on XeMe’ 


This completes the variable selection portion of the 


algorithm, 
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Vil L 


The basic branching method, along with the’ selection 
procedure for which varfable to branch upon are illustrated 
in the following mixed-integer linear programming problem. 
This ts the same problem used in (1] to illustrate the Land 
and Doig procedure. 


Min kx, + 5x 


1 2 
when 5X) + XX) 7 = 2 
x, + 4x, ~ X, = 5 (4) 
3X, + 2X, “Xe 7 


oF > 0, and integer, i =1,...,5. 


The continuous solution produces the canonical form: 


2 t 


+ -X 


= = mi 
1 "10% 7105 “10 


_2 i 
2 aig a WD 


a ae as a 
7 ies». SO 
Zz + 10.4 + 10°5 a 10 ra) 


Since the solution has fractions for at least one of the 
integer restricted vartables, we use the selection procedure 
to tell which variable to branch upon, and which way to 


branch. 
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The computations to select the branching vartlable, and 
the ditirectlon to branch are Included for this’ first 


branching: 


} = ] & S50 89m 1 set ="0.2. 
J = 1 
c¢. = 0, thus no calculation. 
J 
J = 2 
ww Day = 
DL(2) 02 5.59 
XDL =e 5 eo 
Jo= 3 
1 ee I 
= ome = E752 
DU(3) 0k 2 
XDU » 2.175, 

XLOL = 3.5 x 0.8 = 2.8: IML = l, 
XLDU = 2.75 x 0.2 = 0.553; IMU = 1. 
| 2,2 f =@0 0), “ - Fee 0.2. 

J=l 
c= 0, thus no calculations. 
J 
J = 2 
puc2) = Ge! = 2,33, 
0.3 
XDU = 2. 3 oe 
J = 3 
pi¢2) = 44 = 11,0, 
0.1 
XDL = 11.0. 


ADL x fos 11.0 x 0.8 = 388 > XiIepie, 
XDU x (l-f) = 2.33 x 0.2 = 0.467 < XLODU. 
XLDL = 8.8; IML = 2. 


XLDU = 0,55; IMU = 1, 
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>» XLDU. 


f as f 2°092>° 1 - 
J = 1 
C, = 0, thus no calculations. 
y= 2 
iO t 
DL(3) = —— = 2,33, 
0.3 
XDL =" 22535" 
y =3 3 
we) - = = r.08 
erat 
X DU = 1.0. 
MOE x°F = 2.53" x Cl? = 0.467 <-XLoOL. 
XDU x (1-f) = 1.0 x 0.8 = 0.8 
XLDL = 8,8; IML = 2, 
XLDU = 0.8; IMU = 3. 


meow = 8.8 > OFS = XLDU, 


therefore: BSIGN = - 1, 


IM = IML = 2. 


Thus we fmpose the lower bound on 


problem (5) with this added 


form for node 1: 


Xx_> 


» 1, and solve the 


constraint to get the canonical 


wha 


vil co 


5 (6) 


with (24X14 Xq 0X3 Xo Xo) = (35/3, 5/321) 5 3 e).. 
Since some of the integer-restricted variables are stil] 


fractional, we again use the selection procedure to get the 
new bound x, > 2, and solve the new problem taken from (6) 


with the added constraint, xy 22, to get the canonical form 


for node 2: 


Xe - 2x, = 3X4 z= 7 


a - 4x, - Xy ze- 5 
Xs - xg ~ 3X4 z= 2 (7) 
~Zeet OX» + 4x, 2 0, 


where (Zp Xo XoeXaeXyoXe ) = (13,2,1,5,1,1), which is integer 
for all of the integer’ restricted variables, thus, this 
solution becomes the solution vector, SOZZ,Xyreees KX). Next 
the bound A 2 2, ts removed, and the simplex procedure is 
applied to backtrack to the previous node to obtain (6) 
again. Now we impose’ the bound x, £ 1, and solve (6) with 


this new constraint to get the canonical form for node 3: 


Xi + 9X, 7 2Xe = 9 
xX +} oe - 1, = nth 
2 2 1 2 5 2 


x - l. - = (8) 


5 ae 1 a2 
z 2% + 9X 14, 


where (z,x_,X,,X_,X,,X_) = (14,1,2,3,4,0). Since z > 13, we 


lo 2 3 6S 
again backtrack to a node where there is only one branch, in 


this case to node 0, removing the bounds x4 1, and x2 1, 
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and solving by the simplex procedure to again get’ the 
canonical form (5). Then a new constraint of -_ 0, is 
added to (6) and the resulting problem is solved to get the 


canonical form for node &: 


+ = = 8 
x. 10x, 3X 
x + &x - xX = 5 
1 2 ly 
x +lix - x =z 13 (9) 
3 2 ky 
-Z wlilx + 4&x z- 20, 
2 4 
where (z,x ,x_,X. ,X ,x_) == (20,5,0,13,0,8), which again 
1 2 3 h  S 


has z > 13. Since there are no nodes with just one branch 
which can be backtracked to, the problem is solved and the 
solution is: GZzzZ,X% 7X. ,m> xe 7x.) = (13,2,1,5,1,1). Tie 
ice 2 Ss 5S 

procedure is exhibited in the tree shown in figure 1. 

This procedure needed only two nodes to obtain. the 
solution, and four nodes to assure the solution, while this 
problem (4), as shown in [2] requires five nodes. Thus the 


advantages of this procedure over the Land and Doig 


procedure is evident even in this small problem. 
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Node 0 


z= 11.2 
Xo & 0 5 ue ga 

Node & Node 1 

z = 20 z = 11.67 

7 dia) Sass A 
Xs Kk 1 xX 2 ut 
Node 3 
z=ikh 
2 pute 





Node 0: (112/10,18/10,8/10,42/10,0,0) 
Node 1: (35/3,5/3,1,4,2/3,0) 

Nodew2: (33, 2,1,57ie1) 

Node 3: (14,1,2,3,4,0) 

Node &: (20,5,0,13,0,8) 


Figure 1 


Had we just used some arbitrary rule for selecting which 
non-Integer basic varlable to branch upon, the number of 
nodes may have been much greater. The following tree, In 
Figure 2, fllustrates the arbitrary rule of selecting the 


first non-Integer vartable In the basis, x and then 


re 


branching above this varlable by setting the bound, 


x > {x + 1, 
als P 


Node 0 


Node 9 


N oN 
vu 
N = 
N & 


Node 8 


infeasible 





Node 7 


13 
(integer Sol) 
22 = 2 2 13 





Node 5 


z= 17 


(Integer Sol 
2Z = 22:17 





Figure 2 
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Lil 0 10 ULTS 


The algorithm has been programmed in Fortran IV, and run 
on an IBM 360/67. The complete program is” given in the 
appendix. The computer core needed to run a problem is 62K 
plus 36n + 32m + 4m x (m+ n) bytes, where nis the number 
of variables, andm ts the number of constraints in. the 
particular problem when in the form as (1). The’ program 
utilizes a simplex procedure in real variables with a 
round-off correction which requires that al? problems using 
this program have all integer coefficients when in the form 
as (1). 

The results in Figure 3 are for problems in [7], and the 
data for the other programs are from [7] and [8]. 

The results show that this branch and bound algorithm is 
faster than most of the general algorithms for integer or 


mixed-integer linear programs. 
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Comparison of Computational Experience on Haldi and 1 BM 





WO OnNOAMS WHF 


(a) 
(b) 
(c) 
(d) 
(e) 
(f) 


(i) 
(j) 


Source 
ee Ch) 


LIPl | | PM2 | | PM3 BBA BVA 

(a) (b) (c) (d) (e) 
[ter | Iter |Iter | Iter | Time | Time(j) 
(f) (sec)| (sec) 











OumMUUI NWI MN 





NO ND BO ROR 
Ke eS eS OOF roan rf fe & 


Values are from (7). 
Values are from ee 


Values are from {7 
The algorithm presented in this paper. 

The algorithm discussed in [8], times are from [8]. 
Iterations are the total pivot steps required to 

solve the problem. 

The number of constraints. 

The number of variables, excluding slacks. Each prob- 
lem has as many slacks as it does constraint equations. 
An “"F" means that the problem was unsolved after 3000+ 
iterations. 

Times to solve problem on IBM 360 series computers. 


Figure 3 
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ASSEMBLY LANGUAGE TIMING ROUTINE 
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